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Abstract 

The formulation of QCD which contains no divergences and no renormalization pro- 
cedure is presented. It contains both perturbative and non-perturbative phenomena. It 
is shown that, due to its asymptotically free nature, the theory is not defined uniquely. 
The chiral symmetry breaking and the nature of the octet of pseudo-scalar particles as 
quasi-Goldstone states are analysed in the theory with massless and massive quarks. 
The £7(1) problem is discussed. 
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1 Introduction 



In this paper I show how to formulate an asymptotically free theory in such a way that it 
includes perturbative and non-perturbative phenomena simultaneously. 

The idea is the following. Contrary to an infrared free theory, in an asymptotically free 
theory the divergences prevent us from writing even perturbation expansions in a unique, 
well-defined way. We can, however, make use of the fact that divergences in the theory enter 
only the Green's functions and the vertices. On the other hand, knowing the Green's func- 
tions and the vertices we can express all the other amplitudes through them perturbatively 
in a unique way. Thus, we have to formulate equations for Green's functions and vertices 
in a form which does not contain any divergences. If this is done, the solutions of these 
equations will contain both perturbative and non-perturbative phenomena. 

To avoid technical complications, in the first section of the paper we will derive these 
equations in an Abelian theory in which usual perturbation theory is also well defined. 
In the second section we generalize the equations for a non-Abelian theory and discuss 
the connection between perturbative and non-perturbative effects. The equations have an 
integro-differential structure in which the asymptotic behaviour of the Green's functions is 
defined by the boundary conditions. The main conclusion in this section is that in the region 
of large momenta the Green's functions of quarks and gluons contain additional parameters 
compared to normal perturbation theory which can be associated with different types of 
"condensates". These non-perturbative parameters can be defined by solving the equations 
and finding non-singular solutions in the infrared region. A priori it is not clear what type 
of additional conditions have to be imposed on this system of equations in order to fix a 
non-singular solution. It can be, for example, the conservation of the axial current or of 
other currents which is formally satisfied from the point of view of the Lagrangian but is not 
ensured because of the divergences. 

In the third section we will show that these equations make it possible, for the first time, to 
analyse the problem of spontaneous symmetry breaking in an asymptotically free theory and 
to find approximately the value of the critical coupling at which symmetry breaking occurs. 
Also, they allow us to answer the fundamental question for asymptotically free theories, 
namely: how the bound states — the hadrons — have to be treated in these theories and 
how these bound states influence the equations for the Green's functions. 

The analysis leads to the following conclusion. The conditions for axial current conser- 
vation of flavour non-singlet currents (in the limit of zero bare quark masses) require that 
eight Goldstone bosons (the pseudo-scalar octet) have to be regarded as elementary objects 
with couplings defined by Ward identities. This is so in spite of the fact that the couplings 
of these states to fermions decrease at large fermion virtualities. 

The same analysis provides a new possibility for the solution of the U(l) problem. In this 
solution the flavour singlet pseudoscalar rf is a normal bound state of qq without a point-like 
structure. The mass of this bound state is different from zero and can be calculated in the 
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limit of massless quarks. For massive quarks the pseudoscalar octet becomes massive. Their 
masses, however, are not calculable in terms of bare quark masses because of logarithmic 
divergences and have to be regarded as unknown parameters which in the real case of confined 
quarks are defined by the self-consistence of the solution of the infrared problem. These states 
have an essential influence on the equations for the Green's functions, which, as it will be 
shown in the next paper, can be used constructively in solving the confinement problem if the 
effective coupling in the infrared region is not too large. In this case the integro-differential 
equations can be reduced to a system of non-linear differential equations and the theory 
looks like a theory of particle propagation in self-consistent fields defined by the Green's 
functions themselves (as it is the case in Landau's Fermi liquid theory). The self-consistent 
fields are fields of gluons and 7r-mesons. 



2 Equations for Green's functions in QED 

In QED we have two Green's functions: that of the photon D^ik) and of the electron 
G(q), and a vertex function T^k, q). The photon Green's function is defined by the vacuum 
polarization operator H^^k) which can be expressed symbolically by a sum of Feynman 
diagrams 

IW*0 = e^7 M <0 7, + 7„ CD> 7, + • • •} (1) 

In order to obtain series not containing divergences we will try to consider the derivatives of 
II M i,(/c) as functions of momenta. Due to current conservation, we have 

Tl llu (k) = (6 llu k 2 -k ll k u )Tl(k 2 ) (2) 

where U(k 2 ) contains only logarithmic divergences. 

Differentiating @ twice we will have 

d 2 U, v {k) = 6n(fc 2 )^ + (s^ - bj^J (<% + 6)d^U(k 2 ); (3) 

we here use 

oq^ 

The second term in ([3]) does not contain divergences, the first one does. In order to obtain 
a finite expression, consider 

d ll d a U. av {k) = -3^n(fc 2 ) - 3^<9 ? n(fc 2 ). (4) 
Because of this, the quantity 

d 2 U, u (k) + 2d,d a U au {k) = -6^d^U(k 2 ) + U u - ^p-) (<% + 6)^II(fc 2 ) (5) 
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does not contain divergences. 

In Feynman gauge D^ v {k 2 ) is equal 
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1 e 2 (A^ 



-6„ 



k 2 (l-U(k 2 )) k 2 e 2 {k 2 ) ^ v 
Any diagram contains only the product e^D^ik) and therefore 

e 2 D^(k) = ^e 2 (k 2 )5, v 



(6) 



d^U(k 2 



e 2 



(7) 



In first order this means 



1 k^k v 




3 k 



lv la 




the integrals on the right-hand-side are convergent. This fact has to hold in any order. Hence, 
we can include in (|]) the exact electron Green's functions and the exact vertex functions 
and add all the corresponding diagrams which were not included. As a result, we can write 



kudu r, 

e z 



1 h h 



G / \G 





(9) 



where J2 P denotes the sum over the permutations of the indices similarly to @; a quantity 
|^ corresponds to each photon line. Eq. (|9|) is an equation for e 2 (in the form of series in e 2 ) 
provided that G(k) and T^(k,q) are known. 

In order to obtain equations for electron Green's functions and vertices we have to re- 
member that these functions can change rapidly even if e 2 is small because they can have 
infrared singularities of the type e 2 In 9 ~™ or ultraviolet singularities of the type (j^j ■ It 
is proven to be possible to arrange the differentiation in such a way that there will be an 
expansion only in e 2 . 

To write an equation for the fermion Green's function not containing any divergences, 
we have to differentiate twice the self-energy of the fermion or its inverse Green's function. 
Let us consider d fl ^ l G~ 1 (q); the diagrammatic expression for G _1 will be the following. The 
simplest diagram is 
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Diagrams of the next order are 



--- k k 



It can be easily shown that 

d 2 — — = -ATi 2 i5\k) (10) 
k 2 + «e 

Using this equality, we have 

_^ e 2 



<9 2 ^7 ^7" = -^o7m G i7 m > 



,2 — 



Restricting ourselves to d 2 -^, we can write 

dVi = -^oT^oTau 

5V 2 = -go'j^God^ - god^Goin - go^G^ - gi^G ^ 
and, consequently, 

^G- 1 = gd li G- 1 Gd lt G- 1 , G- 1 = Gq 1 + G^ 1 + G^ 1 . (11) 

The term 5 2 has a contribution to , " 7 v ~~ \ not containing overlapping divergences and 
is of the form 

_2fc CT _iK 

7a 7a 

_ 2fc ff _ 2fc^ 

+ + (12) 



7a 
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For the calculation of higher order diagrams it is convenient to adopt the following principle. 
Beginning from the first point of interaction we shall relate the external momentum to the 
photon line: 

(13) 




q 

The next photon interaction can be expressed as 





(14) 



now we relate the external momentum to the positron line. If the positron is emitting a pho- 
ton, we keep sending the external momentum along the positron line up to its annihilation. 
As a result, we get a tree structure 



(15) 



in which the initial photon line can end only at the final electron. It can be shown that 
every line in the diagram will be passed only once if the photons included in the fermion 
self energy are not taken into account. This means that in this approach exact electron 
Green's functions have to be used with bare vertices since the idea is basically the exclusion 
of overlapping divergences. Taking the second derivative of one of the photon propagators 
and restricting ourselves to the contribution Air 2 5 i {k)g{0) we obtain both on the left-hand 
side and the right-hand side photon emission amplitudes with zero momentum. Due to gauge 
invariance, however, a zero momentum photon cannot change the state of the system (it is 
emitted from the external lines). The photon emission amplitude equals T M (g, 0) = d ll G~ 1 (q) 
[within this convention the bare vertex is — cy^]. The final contribution of the differentiation 
has the structure (|i~TD. The above statement, which essentially means that the emission of 
a zero momentum photon is determined by the total charge, can be proven by the Ward 
identity. 

Thus, we have 

d 2 G- 1 = g^d.G^Gd.G- 1 + -+Jf*2£^y+- (16) 

The remaining diagrams contain first derivatives of photon and positron lines and second 
derivatives of positron lines. All these diagrams can be expressed in terms of the exact T M , 
G and D functions. They do not contain divergences except those graphs which correspond 
to the photon self energy. 

The first term in ([16]) has the structure 

±g(0)M w (q,k = 0) 
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where M vu (q,k = 0) is a quantity close to the Compton scattering amplitude in the sense 
that they would be equal if we differentiated all fermion propagators including those inside 
the electron Green's function. (The factor | is due to the fact that the Compton amplitude 
contains the sum of diagrams with momenta k and —k). The Compton scattering amplitude 
M uu (q, k = 0) satisfies the Ward identity 



M uv {q) = G- 1 d 2 G G- 1 = 2<9 At CT 1 G d^G' 1 - d 2 G~ 1 



which has the simple diagrammatical meaning 



fc = 




/ k'=0 



+ 



(17) 
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The second term on the right-hand side of ( fL7]) represents the contribution of photons which, 
if we carried out the mentioned differentiation, would correspond to photons inside the 
Green's function and which are not present in M w . Hence, 

M yv (q^) = 2d tx G- 1 Gd^G- 1 . 
All the diagrams for d 2 G~ 1 are of the form 




and contain the photon lines with the exact Green's functions (H) which equal 47r 2 g(fc)p-. 
By differentiating the photon lines we have calculated the contribution of d 2 -^. So, there 
remains 

f d*k M vv {g,k) ^ g(k)-g(0) 

I -~J4^H 2 9 P ' (19) 

where we have introduced g(0) in order to avoid a contribution from 5 4 (k). This expression 
contains logarithmic corrections coming from the ultraviolet region (k > q). Replacing 
g(k) — g(0) by g(q) — g(0) + g(k) — g(q) and performing an integration by parts, we get 



1 



Muu(q, 0) 



d*k g{k)-g{q) d2 M uu (q,k) 



4ir 2 i 



k 2 



(20) 



The first term in ( p0| ) can be explicitly calculated while the second one does not contain any 
logarithms because of the presence of the difference g(k) — g{q). We can rewrite (|20D in the 
form 

X = (g(q) - gmd.G-'Gd.G- 1 + 6 i; 
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consequently, d 2 G 1 can be expressed as 

&G~\q) = g(q)d,G- 1 Gd fM G- 1 + S, + 5 2 (21) 

where 5± is defined by [the second line of] fl20D and S 2 contains first order derivatives of the 
photon and positron lines and second order derivatives of the positron lines, as it has been 

2 2 / \ T 

explained above. The first term contains all singularities of the types a In g ~^ and ( — ) . 



3 Equations for the vertex function and for the ampli- 
tudes of interaction with the external field 



To obtain the equation for the vertex function, let us express ^(p, q) as a set of diagrams 
containing exact Green's functions: 



T^p,q) 





4+1 + 




+ ... 



(22) 



We shall relate the external momenta to the lines in the diagrams in the same way as we 
did when we derived the equation for the Green's function. Calculating the second order 
derivative in q we obtain 



%r»(p,q) = -g(0)MH(p i q,k) 



k=0 



+ ... 



(23) 



where M£ u (p, q, k)\ k=0 is the contribution to the amplitude of the process 




which corresponds to the emission of photons with momenta k\,k 2 = from the external 
legs: 



AC(p,g) 



(24) 



Inserting into (|23"D and replacing g(0) by g(q) we get 

d 2 T^p,q) = g(q){A v (q 2 )d v r»(p,q) + d v r»(p,q)A u ( qi ) 
- A u (q 2 )T»(p,q) A u ( qi )} + d 2 f»(p,q) . 



d v G~ x 



(25) 
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Here we have introduced the notations 



gi,2 = g±|, A,(q) = d,G- 1 (q)G(q) and A^(q) = G(q)d fl G^ 1 (q). (26) 

The correction terms <9 2 T M are defined as a set of diagrams with exact Green's functions 
and vertices. They contain first order derivatives of both the photon lines and the positron 
lines, second order derivatives of the positron lines and, in the same way as in the case of 
the Green's function, corrections due to the replacement of g(0) by g(q). 

The same equation is valid for the interaction amplitude of fermions with the external 
field provided this interaction does not depend on the relative momentum q of the fermions. 

The equations for the interaction with external fields are essential. Indeed, if these 
equations have solutions that decrease at large virtualities of the fermions — i.e. solutions 
which do not require driving terms — this means the existence of bound states. 



For the sake of convenience we rewrite the equation in the form 

d 2 <j)( P , q) = g{q) {A„(q 2 )d u <P(p, q) + d u (j){p, q)A l/ (q 1 ) - A u (q 2 )0(p, q)A v ( qi )} , (27) 

which will be understood as the equation for the bound state with a spin which is defined 
by the invariant structure of the matrix 0. 

It is important to note that the accuracy of the equations for the vertex T^(p, q) and for 
the bound states differs from the accuracy of the equation for the Green's function. The 
functions T^(p,q) and 4>{p,q) depend, among others, on the ratios p 2 /qf, p 2 /q 2 - If these 
parameters become large, then T M and <p(p, q) contain, in general, the so-called Sudakov 
logarithms which are not included in the equations (EH) and (27). Hence, the equations are 
valid only if 

a In In ?- < 1 . (28) 



4 Equations for Green's functions in QCD 

QCD is the theory of interacting quarks and gluons. The description of quarks is more or 
less the same as in QED. Gluons, however, are very different. Even the fact that being a 
spin 1 particle the gluon has to have a multi-component wave function is not seen explicitly. 
In the usual approach it is described by a Green's function D^ u = (A^(x) , A u (y)) which in 
momentum space (in Feynman gauge) always can be written in the form 

= ^C(k?). (29) 

It contains only one unknown function. All spin properties of the gluon are included in the 
momentum dependence of its interaction. In order to introduce multi-component Green's 
functions, we have to formulate the theory in a form in which the interaction is momentum 
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independent: we have to replace the usual description of the gluons and their interactions by a 
Duffin-Kemmer type formulation. In the framework of this description the gluon Lagrangian 
is 

C[x) = -{dpA, - d u A^ + g[A u , A^F^ + (30) 

[where the trace over colour indices is implied]. The potential A^ and the field strength 
are independent quantities. The interaction is momentum independent and equals 
[Afj,, A^F^y. In this formulation we have three independent Green's functions 

{A^A V ), (F^Ap), (F^,F pa ). (31) 

In order to fix the gauge in a covariant way, it is, of course, necessary to introduce ghosts 
by adding a gauge-fixing term 

V = \{d,A,f. (32) 

The three Green's functions can be combined into one by introducing the ten-component 
state [implying p < a] 

*=(/;)■ ( 33 ) 

We use the parameter m, which has the dimension of mass, to convert the lower compo- 
nent of \I/ into the same dimension as the upper component. The equations for the fields 
corresponding to the Lagrangian (|30|), (|32|) are 



V„F^ + Cd»d v A v = 0, 

d^A v - d u A p + g[A^ A v \ - = 0. (34) 
In terms of the state \l/ we have 

^ (idp + gA p ) m- mi ^- l±C(p 2 - p 2 )q = 0; (35) 

m 

p = dfj,, p 2 = -P^Pudf.du , p 2 = -d^. 
Here (3^ are Duffin-Kemmer matrices satisfying the commutation relation 

PiPkPi + PiPkPi = 5 ik pi + 5 ki pi . (36) 
They are connecting the [four] upper ty" and [six] lower \l/ p(T components of the field \& in 














[Their non-zero components, the 4x6 matrices S^} have a simple representation 

(StiYpa = i{$W$v<r - <W<^p)- (37) 

The quantities j± are projection operators onto the upper and lower components of \I/ . 
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If we want to preserve the dimension of the (A^, A v ) component of the free Green's 
function, D has to satisfy the equation 



P_ 

rn 



7- 2-({P -P ) 



D, 



1 



— 2' 



The solution of this equation is 



Do 



P 



m m z 



(38) 



(39) 



In Feynman gauge ( = 1 and C\ = 1, in Landau gauge C\ = The three terms in ( p9[) 
correspond to the three independent Green's functions (pTj) ■ 

The vertex for the interaction of three gluons with momenta Pi,P2,P3, colours a, b, c and 
Duffin-Kemmer indices a, (3, 7 has the structure 



a,a,pi ;6,/3,p2;c,7,p3 



P2 



P2 



Pi 



+ 



P3 



Pi 



+ 



P3 



Pi 




P2 
P3 



(40) 



+ 



/./, 



v pa pa v 

The coupling constant remains in our notation g. 

Let us consider in this approach the properties of the exact Green's function 
IT 1 = k - 7- - 7+ C{k 2 - k 2 ) - £ - 5(Ck 2 - k 2 ), 

P 



k 



m 



where E is defined diagrammatically: 





(41) 



(42) 



It contains four matrix elements £++, £+-, S It can be also represented in the 

form 

S = fcS 1 +7_S 2 + 7 + Ps 3 - (43) 

The factor k in the third term is necessary to preserve the current conservation; in the first 
order the ghost contributes only to S 3 . 

We added the term 5( in fl4l"l) to be able to fix the gauge for the exact Green's function. 
Instead of fl4"T| ) we can write 



zr 1 = z^k - z 2 -V - z^k 2 1+ - c(P - k 2 ). 



(44) 
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The Green's function then will be equal 

1 



D 



(z r 2 - V 



{z- L k + Z~ l C 1+ + Z 2 Z~ 2 k A ^} + 



Z 2 7- 



where 



C= 1 + 



k 2 



z 2 z x 



zi 



In Feynman gauge 



-2 



(45) 



(46) 



(47) 



In order to understand the meaning of the denominator in (45), let us consider the renor- 
malization properties of simple diagrams, for example S + _ + 




(4£ 



where g m is the effective bare coupling; [we expect, see ( |56"D below] 

r+-+ ~ Zf 1 . 

If this is true, we have in Feynman gauge 



+ E. 



Z^Z 2 l 



(49) 



(50) 



However, p A1 <9 M XL + [in the manuscript, S ] has to be proportional to a.Z x . This means 

that we have to expect 

W _«(p) x (51) 



i.e. 



Z 2 Z X 



-2 



^z r 2 z. 



This is our definition of a(p). It has all known properties of the renormalized coupling a. 
As a result, D can be written in the form 



D 



'a(p 2 
a 



] -ZiZ 2 l {z 2 Z^k + C 1+ + (Z 2 Z^k) 2 j-} - 2 + . 



(52) 



In this approach is not a quantity coming from the normalization: it is the bare coupling. 
The theory has to be defined as the limit «o — > 0. In this context the expression (|52|) has 
a very interesting property. In the limit ao the Green's function D contains only Z{ 1 
and Z 2 . According to (51), in this limit Z% = Z± Z 2 . Because of this, when we regard 
the equations for Z-f 1 and Z 2 X in the way we did for the fermionic Green's function in QED, 
ao has to disappear. Consequently, we have equations for Z-f 1 and Z 2 X with a(p 2 ) being 
arbitrary. The equation for Z-f 1 will not help since, due to the equality Zjf 1 = Z{ 2 Z 2 it 
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has to be an identity. We will see that an equation for a appears when we will consider the 
correction of the order of tJolq. 



Before formulating the equation for the Green's function, let us see what the Ward 
identity looks like in this formulation. The bare vertex is 

rj = M, 

[with / the colour matrix (the structure constant)]. Consider the relation between (p 2 — 
Pi)^ >l {P2,Pi) and the Green's function: 



Pi*- Ph 



P2 - Pi = m 



1- - {k - 7-) 



= m[D \k 2 ) - D^(h)} + m 7+ Co[(A; 2 2 - k\) - (k\ - k\)\. 
Hence, at p — > [hereafter, = d/dp^] 

T°\ p= o = f Imd^D, 1 + m 7+ Co^(P - k 2 



(53) 



(54) 



which is, of course, the usual complication due to the Slavnov- Taylor Ward identity. But in 
this formalism the additional term is equal 



-mdft Cc 



d_ 

Wo 



D. 



-i 



Inserting this vertex into an arbitrary diagram we obtain [for the full vertex] the relation 

d 



fm 



duD- 1 - CoTz-d.D- 1 



According to (0), (|47|), in Feynman gauge we can write [d^ = ( d/d( } 

= }md, Uz^-d c z^)k - 7 _(z 2 - 1 -a c z 2 - 1 ) - 1+ k 2 (z^ 2 z 2 -d c (z^ 2 z 2 )) 



(55) 



(56) 



The quantities d^Z{ x , d^Z 2 x must be calculated from the equations for Z{ 1 and Z 2 X . But 
if the equations are formulated in terms of the exact Green's functions, the dependence on £ 
enters only through these Green's functions D which, according to fl4*5"D , ([46|), include £ only 
through the quantity C defined in fl46|) . In the limit a — ► 0, however, C does not depend 
on (. Consequently, c^Zf 1 and d^Z 2 are equal to zero and we have the following simple 
relation for the vertex at zero momentum: 



Here D~ l is defined by (|56D; it does not contain gauge-fixing terms: 

L)- 1 = Z^ l k - 7_Z 2 ~ 1 - 7+Zf 2 Z 2 k 2 . 



(57) 



(5? 
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Let us first consider the equations for Zf 1 and Z 2 X in the same way as we did for fermions. 
As in the case of QED, we will use the Feynman gauge in order to simplify the one-gluon 
contribution to the equation for the Green's function. 




<X> +-..} (59) 



Each line here contains the exact propagator (p2|) . Differentiating the propagator along the 
upper line, we obtain the contribution corresponding to the second derivative of one line 

d 2 D = -47r 2 z5 4 (p) 7+ ^(0) - ^ (NZ*) + ^d 2 (HZ?) + d 2 ^, (60) 

where 



,2 7 -l«(°) 



Z = Z{Z 2 '-^ (61) 



and we introduced the notation 



Zi Z<vi^ 

D = —N+^ r . 62 

The contribution of this derivative to the equation will have the form 

- 4:7i 2 iZ^M{p,p) +M V (63) 

Here M(p,p) is the gluon-gluon scattering amplitude at zero momentum of one of the gluons. 
The second term Mi is defined diagrammatically: 

M 1 = 



Taking the first derivative of two different lines, we get the contribution 

M 2 = -<^^>- + .... (64) 

As a result we have 

d 2 (Z^k - Z^j-) = —Z^M(p,p) + Mi + M 2 . (65) 

7T 

We have to remember that on the right-hand side we have to take only the matrix elements 
H+)> (+|-) and (-|-)- Writing 

M(p,p) = r M D T M + M(p,p) (66) 
and using the Ward identity (|57|) , we obtain an equation of the same structure as for fermions: 



0> { z?i - z 2 -S-) = /^^!M ejri L + M±,_ | a ^ + l, (67) 
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L = M + M 1 + M 2 , f 2 = N c = 3. 

If we now make the same trick as before, namely replace Z(0) by Z(p) and redefine L, we 
will have in the limit cto — » 

d 2 Z^[j_, k] = [ 7 _, d^D-^Zyk + ZlZ^ 1+ + Z^V)^- 1 ] + [ 7 _, L'}, (68) 

9 2 Z 2 -S- = 3^ l-d^D- 1 (Zxiz + Z\Z^ 1+ + Z 2 £ 2 7 _) ^D-S- + 7-^7- (69) 

7T At v 7 

[Equation for quark Green's function.] The equation for the fermionic Green's func- 
tion in QCD will differ from ([2lD only by the factor A a A a = | (A a are colour matrices). The 
reason is the following. The derivation of the equation (f2lD was based on the relation be- 
tween the fermionic Green's function and the amplitude of zero momentum photon emission 
by a fermion 

r M (g,0) = 9 M G- 1 (g). (70) 

In the usual formulation this relation is not correct in QCD. The simple relation fl5"T|) for 
the amplitude of the zero momentum gluon emission by a gluon implies, however, that the 
amplitude of a zero momentum gluon emission by a quark has to be equal 

r;( q ,o) = x a (z^z 2 )^G-\ q ). (71) 

Together with (|52| ) it leads to the equation (|21|). 



[Equation for T.] The equation for the colourless vertices remains also the same. 

The equation for a three-gluon vertex, however, will be essentially different. In the same 
way as for the fermionic case, but taking into account the non-commutativity of the gluon 
coupling, we can show that 



vppu 



lap 



satisfies the following equation: 



d 2 t» [q+ P -,p, q -^ 

- A u ( Pl )T»A u (p 3 ) - A u ( Pl )t»\A u (p))^ - {A v {p))^'A u {p 3 ) }; 



V 



(72) 



A u = d„D- l D 



A„ = Dd„D 



-l . 



Pi 



, V 
q+ 2 



p 

P3 = q ~2- 



The right-hand side corresponds to all possible gluon emissions from the external lines: 



+ 



+ 



d v D- 



+ 
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+ 



+ 



+ 



d„D- 



dvD- 1 r dvD- 1 d v b- 



b v d- x r 



Higher order terms can be written in the same diagrammatic way. 

The fourth term in (|72|) and in the corresponding diagrammatic expression equals zero 
at a = since DTD = at a = 0. 

The equation ( j72|) is, indeed, quite different from the equation for a colourless vertex. 
The main difference comes from the fact that on the right-hand side it contains derivatives of 
T not only over q but also over p. Hence, in order to find T, we have to write three different 
equations for second derivatives over three different momenta. 

A similar equation is valid for the quark-gluon vertex. 

Knowing the Green's function and the vertices, all the other amplitudes for the interac- 
tions and quarks and gluons can be written in the usual perturbative way. These amplitudes 
have no divergences and contain inside the gluon Green's function the unknown function 
a(p). To formulate the theory in an unambiguous way, without any references to the cutoff 
and the regularization, we have to find the equation for a(p) as we did in QED, and learn, 
how to write the higher terms more elegantly and constructively. I postpone the investiga- 
tion of this problem to another paper. What I want to discuss now is the main difference 
between an asymptotically free theory and an infrared free theory. 

The equations for Green's functions of quarks and gluons are proven to be second order 
integro-differential equations. To solve them, we need boundary conditions. In an infrared 
theory the boundary conditions are known: they are the conditions for the existence of free 
elementary particles at small momenta. In an asymptotically free theory the interaction is 
small at large momenta, and we expect to have here a perturbative solution. However, in 
the region of large momenta all the equations have two types of solutions: the perturbative 
solution and solutions which decrease as some power of the momenta and therefore contain 
dimensional parameters reflecting the density of different condensates. These parameters 
have to be defined either by the introduction of additional conditions of the type of con- 
servation laws or by the self-consistency of the solutions in the small momentum region (or 
both). The next section of this paper will be devoted, essentially, to the discussion of this 
problem. 
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5 Spontaneous symmetry breaking in asymptotically 
free theories 



In the present section we will consider the equation for the fermion Green's function in QCD 

d 2 G-\q) = g(q)d,G-\q)G(q)d,G-\q) , (73) 



where 



g{q) 



a(q 2 



-X a X a 



4a(q 2 



71 37T 

This equation was discussed extensively [1] in connection with the problem of quark con- 
finement. Here we shall use the equation for the discussion of the spontaneous breaking of 
chiral symmetry in asymptotically free theories. The asymptotic freedom is reflected in this 
equation by the fact that a(q 2 ) decreases when q 2 — > oo. We will assume that in the limit 
q — > a{q 2 ) approaches a finite value. 




Figure 1 

At q 2 — > oo the solution of the equation fl73p has the form 



G-\q) = Z- 1 



(m- q) + — + 



r 



^q 

q 4 



(74) 



q = i»q» ■ 

If a = 0, the quantities Z, m, v\, v 2 are arbitrary parameters. In the case when a(q 2 ) is 
defined by perturbation theory, Z,m, V\, v<i are the following functions: 



Z 



OqJ 



a \ 7 " 
m = m | — 
a J 



Via 



(75) 



The anomalous dimensions jz, 7m and 71^ can easily be found from the equation (|73|). 
Generally speaking, the solution depends on four parameters. In the limit q 2 — > 00 the 
chiral invariant solution can be written as 



G 



-1 



Z- l q I 1 



(76) 



The general solution ( |73D corresponds to massive quarks. In the solution which corresponds 
to spontaneously broken chiral symmetry, m = 0. In this case the mass term decreases 
when q 2 — ► 00; the term v\ is responsible for the violation of the symmetry. 
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Multiplying ( [75| ) by G(q)G 1 (q) we obtain the equation 

d 2 G- 1 (q)=gA ll (q)A ll ( q )G- 1 (q) (77) 

where 

A^q)=d,G-\q)G(q). (78) 

Clearly, it has a structure which corresponds to the equation for particle propagation in the 
self-consistent field gA^A^. It is easily seen that (|73|) is equivalent to the equation for A^ of 
the form 

d^A^ = -PApAp, P = l-g. (79) 
The matrix G~ x is defined by two invariant functions and can be written as[] 



G- l = Z-\q)[m{q)-q)=pe 



-nf 



10) 



where 



Here A^ is 



n = - , q 

q 



A 



- — ~ \ndn<t> - d^n sinh^e n 2 



for a constant (3, or 



it 

p = — exp 
9 



£ ("It - A -L 



(3 \u 



for (3 which is a function of q. Here 



51) 



Inserting (|80|) and (j8l| ) into (|77|), we obtain a set of non-linear equations for p and (p. We 
can linearize the equation for p by writing 



(82) 



(83) 



du du 
u = q U Tr- = quOvii = — , 
oq y <9£ 



£ = \nq. 



As a result, we get for u and (f> the following set of equations: 



il-u + (3 2 



3 sinh — h 



u • 

+ 2—0 — 3 sinh < 



0. 



For a constant /3 the conservation law 



(84) 
(85) 



dt-E = 



1 [p can be treated as a dimensionless quantity since the equation uTjn is homogeneous] 
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is fulfilled; 



E = u 2 - u 2 + (3 2 



3 sinh 



u 2 . 



*6) 



The term ^{u — u) is of the order of g 2 and it can almost always be neglected. 



The asymptotic behaviour 
corresponds to 



J) of the Green's function in the limit q 2 

u — > Cq 2 , q 



111 . 



The chiral-invariant solution corresponds to = m. Close to the value 
in + 0) , at large q 2 we have 

u 



oo {(3 



in i.e. 



a 



1 + 3/? 2 . 
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Hence, fl8"5|) is an oscillator equation with damping if q 2 is growing and with acceleration if 
q 2 is decreasing: 

+ 30= -2^1 + 3/3 2 0. (88) 

This means that the chiral invariant solution $ = in is unstable in an asymptotically free 
theory. 



Let us consider the equation for at negative q 2 values in detail. Due to (PD|), = iip is 
in this case purely imaginary and ([85]) describes the motion of the particle as a function of 
the "time" £ in a periodic field; the damping (or the acceleration) is defined by 



1 














— »- 











Figure 2 



At £ — ► oo (|g 2 | — > oo) the particle is situated at one of the minima of the potential; it 
accelerates as £ decreases. The acceleration rate is defined by the parameters of the solution 
(0) in the limit \q 2 \ — > oo. When q is decreasing we have, generally speaking, two possible 
behaviours for the solution. 

If - remains positive all the time, then the solution goes to infinity and G develops a 
singularity at some Euclidean q. If instead the potential - changes sign, then the solution 
may approach again a minimum. In the latter case it is easy to see that G~ l has a singularity 
at q — 0. 

There is only one possibility to avoid having a singularity in the Euclidean region includ- 
ing q = 0. We have to choose the parameters which determine the acceleration at \q 2 \ — ► oo 
so that the particle approaches the maximum of the potential when q — * (£ — ► — oo). In 
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order to find these exceptional solutions without singularities at q 2 < it is natural to start 
solving the equation at q = by fixing the behaviour of the solution at £ — ► — oo. The 
solution of (184]) ) fl8~5l) corresponding to a maximum (e.g. ip — 0) at q — > is 



£■0 



u = Z- l (Q)q 



(89) 



In the case of such a solution the constant E in (|86f) equals zero, and 



u 
u 



\ 



1 + /3 2 



o • 2^ ^ 

3 sin — 



(90) 



Inserting ( p0[ ) into (pS|) we obtain one non-linear equation. It can be analysed easily for 
arbitrary q 2 values. The solution we are interested in contains essentially one parameter m c 
which can be related to the renormalized fermion mass. The parameter Z^ 1 is irrelevant; 
it defines the normalization of the Green's function at q = and can be chosen as 1. 
The solution of the equation leads to the unambiguous determination of the asymptotic 
parameters of the Green's function ( |74|) by m c and by the parameter Aqcd of the strong 
interaction^ which enter (3{q) = 1 — Cpa{q 2 ) /it. 

As it was mentioned before, to the spontaneous breaking of chiral symmetry corresponds 
an asymptotic behaviour of G _1 in which m = 0. The existence of such a solution requires 
a connection between m c and Aqcd- The renormalized fermion mass as well as the other 
"condensate" parameters are then determined by the strong interaction parameter Aqcd- 

Let us consider this solution in detail. Starting from the point ip = at £ — > — oo, 
the solution will, obviously, reach the minimum of the potential either monotonically (if the 
damping is strong enough) or in an oscillating way. In our case the damping depends on 
the value of (3. If g is small, (3 is close to unity, the damping is strong and the solution has 
a monotonic behaviour. By decreasing (3 the solution may become an oscillating one. For 
obtaining the value of (3 at which the oscillation begins there is no need to solve the equation 
for all values of q. It will be sufficient to investigate the solution in the region where ip 
becomes close to 7r; here ip = n + ip and ip satisfies the equation ([88]). The solution can be 
written in the form 



It oscillates if 



= e - K Ccos [J2- 3p 2 £ + 5 



(3 2 < 



9c = 1 



3 3 V 3 

Oscillations in ip mean that the mass term in (|8( 



V 



< 9 < 1 



1 + 3/3 2 




(91) 



(92) 



m(q) 



i— q , 
2 H ' 



(93) 



2 [Aqcd is the momentum scale where the coupling in Figure 1 becomes of order unity or, more precisely, 
when it reaches the critical value (see below)] 
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starts to oscillate. By solving the equation for bound states we can check that at g > g c bound 



states appear with wave functions behaving like (|91~D which coincides with the behaviour of 
the solution of the Dirac equation in the field of a point-like static charge Ze when Z > 137. 
The simplest example for such bound states are Goldstone states the wave functions of which, 
as we shall see in the next section, are proportional to m(q). 



We have found the oscillations and the critical value g c — 1 — y~ using the assumption 
that g = const. In reality g depends on q (see Fig. 1) and the situation is somewhat more 
complicated. It reminds the case of the equation for a critical charge of finite radius. In the 
region of small q values where g(q) is close to a constant g(q) « g(0) we can consider ip as 
an independent variable and q 2 as a function of ip. It can be shown that for g(0) satisfying 



the condition 
plane (Fig. 3) 



there are two regions < ip < tt — if) c and tt + if) c < ip < 2ir in the q , ip 




Figure 3 

where the solution if)(q) is monotonic, and there is a region tt — ip c < ip < tt + ip c where the 
solution oscillates. The value of ip c is determined by the equality 

2 A (i 




Sin 2 = If ^ ^ : (9 U 

it exists if (3 2 < |. Considering [3 as a function of q 2 in ( p4[) and taking into account that 
f3 2 — > 1 at \q 2 \ — > oo, we obtain a region bound by the dashed curve and the vertical axis 
where the solution oscillates. There are no oscillations if \q 2 \ > Aq CD (/? 2 (Aq CD ) = |). Due 
to ( [741) and (|75|) , we can write in the region \q 2 \ 3> Aq CD 

2 g Va / q \a J 



If /3 > |, the solution which equals iip = q/m c at q —> 0, transforms monotonically into ( p5|) 
with mo 7^ 0. If f3 2 < |, ^(Aq CD ) and V'(Aqcd) start to oscillate as functions of m c , and m 
can turn into zero at a certain m c value. This means that we have a solution corresponding 
to broken chiral symmetry. 

If mo = 0, there exist also a large number of solutions depending on the parameters v\, 
z/ 2 . With the same sign of v\ we can have different solutions corresponding to the curves 
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I and II which in the limit q — > reach tp = and ip = 27T, respectively. The solutions 
V>(0) = 2?t correspond to smaller values u±, m c . 

Let us consider the solutions of types I and II in detail for complex q — this is justified 
since G~ x has to satisfy the requirements of analyticity and unitarity. We will show that 
both solutions have singularities at real positive q 2 values. The solutions are chosen in such 
a way that they are regular as q — > and have no singularities at q 2 < 0. Due to the 
analyticity of the equations, the behaviour of the solution for q 2 > can be found by solving 
the same equations (|85|), (|90| ) with the same boundary conditions at q — > 0. 

If (3 is fixed, the equations flgip - ( IgCT ) can be rewritten in a simpler form. Denote 



u 



p(</>). 



u 



Then 




(96) 




(97) 



with the boundary condition p = 1 



at = for a I-type solution. 



[ Equation (|9j]) solves flSBJ) for E 



0. Equation ( ^6] ) follows from ( P7| ) and 




which is equivalent to eq. ( 84j) for /3 = 0. ] 



The phase diagram corresponding to the equation ([96|) is shown in Fig. 4 where the solid 
line represents the solution of the equation 




The function p = p((f>) is shown by the dashed line in Fig. 4. 



/ / \ \ 

/ / \ \ 




Figure 4 
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The ^-dependence of p at <fi — > oo is different for ft > | and ft < \. In both cases 
approaches infinity at finite £ values. We have at ft > |, — > oo, 

§|=/%>, P = ~ e»; C>0, (99) 

and at ft < |, — > oo, 

p = -2/3 2 C ie 2, 



I 



° l 2 V 1 - 4/? 2 
The equation (|97|) enables us to find £ as a function of 0. 

ft 



s ^ 2 A 



p 2 (0) + 3/3 2 sinh 2 §- 1 



(100) 



(101) 



The integral in ( |101| ) converges in both cases (/3 > | and ft < |); because of this, goes to 
infinity at a finite £ = (g — > m*). Near q = m* at ft > ^ (g < 2) we have 



tt = itoWl 



e 2 



rrr 



m* 



and, consequently, 



G-\q)=Z^ 1 



+ g / 1 \ g — g 



+ 



C 



If /?<§ (<?>§), 



It = Uq I 1 



2/3 2 



= d 1 



77T 



and thus 



g \ 2/3+1 g _|_ £ 



g ^ 2/3 1 1 g — g 



rrr 



C{ 2 



(102) 



(103) 



(104) 



(105) 



It can be shown that the relation between the position of the singularity at q — m* of the 
Green's function and the quantity m c can be written as 



rrr 
In — = 
m c Jo 



ft 



2wV(0) + 3/3 2 sinh 2 f - 1 sinh20 



(106) 



The formulae (|103|) and ( |105| ) define the behaviour of G 1 (q) near the singularity for a 
solution of type I. 

In order to obtain the behaviour of G~ l near the singularity for a solution of type II it is 
sufficient to notice that at q — > such a solution has the form 



2m - X. 

m' 



(107) 
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The replacement of by + 27r* changes only the sign of Replacing q/m c by —q/m' c 
and solving (|96| ) and (0) for g > we will find to be negative, and near the singularity 
will go to — oo. This means the following. If the first term G7 1 (q) of the expression 



G-'( S )=e;'( 9 )i(i + f) + e:'to)i(i-f) 

equals zero for the solution I at q — m*, then Gz 1 is zero for the solution II at q = m 
If q 2 > m* 2 , both solutions become complex: 



(108) 



> 

-2 In 



_q_ 

m* 


A 


+ 


VK 
J 


_Q_ 

m* 




+ 


2m 



'4 



(109) 



The trajectories of 0(g) at < q < oo in the complex plane are shown in Fig. 5 for the 
solutions I, II [for f3 > |]. 



ir 



ii 



2Z7T 



ITT I 



Figure 5 

The solution of the type I has a remarkable feature: Im > 7r for any q > m* values. Hence, 
taking = 4>i + i<p2, the imaginary part of m(q) in (pOj) 



Imm(g) = glm 



-g sin 02 



cosh |(0i + Z0 2 ) 
sinh |(0i + z0 2 ) 2| sinh |(0i + z0 2 ) 



(110) 



turns out to be positive. At the same time Imm(g) for the solution of the type II is an 
oscillating function; this can lead to a contradiction with the unitarity condition for the 
Green's function. 

For complex q values G _1 (g) has no singularities. If we move in the complex plane along 
an arbitrary ray [from the origin to infinity], the trajectory of 0(g) doesn't approach infinity 
(curves I' and II'). 

The main result of this section is the following. In the framework of the equation for 



the fermion Green's function (73) there exist solutions corresponding to broken symmetry 
(q. 

behave at q 2 — > oo as 



provided g(q) has an asymptotically free behaviour and g(0) > 1 — y|. These solutions 



G- L (q) 



Q \ Q J 
The expression ( |1 1 1[ ) has a mass term decreasing at infinity. 



(Ill) 
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6 Axial current conservation and Goldstone states 



If a fermion Green's function corresponds to symmetry breaking, it is natural to expect the 
existence of Goldstone-type bound states. This expectation is connected with the belief 
that if m = 0, the axial current has to be conserved. This is, however, not necessarily 
true because due to divergences in the theory a leakage of the current is possible into the 
region of the ultraviolet cutoff. A typical example for this phenomenon is the anomaly. But 
even in a non-anomalous case it is not obvious whether the current conservation is implied 
by the equation for the Green's function or it is imposed as a condition on the solution of 
the equation. In order to clarify this, let us consider the equation for the bound state 
supposing that it is a pseudoscalar. 

<9 2 0(p, q) = g(q) [A^q 2 )d^(p, q) + drffa ?)i M (?i) - A^q^p, q)A,( qi )} . (112) 

This equation has to have a solution decreasing for large q 2 at p 2 = 0. It is easy to see 
that indeed there exists such a solution. At p = the equation for is an equation for the 
variation of a function which satisfies the equation for G^ 1 . If this variation is taken in the 
form = C{7 5 , G -1 }, obviously satisfies the equation and decreases at \q 2 \ — > oo as 

- C 2 -4 lb . (113) 

This, however, does not mean that we have particles with p 2 = 0. Indeed, the equation 
( p.!2[ ) has a solution decreasing at \q 2 \ — > oo for any p 2 values. The reason for this is that 
the equation is highly degenerate. It has a solution of the form 

(f ) = O l G-\q 1 )+G-\q 2 )0 2l (114) 

where 0\ and 2 are any combinations of Dirac matrices with coefficients not depending on 
q. This can be checked by substituting directly (|114j) into ( |112|) and using the equation for 



G^ 1 . The reason for this degeneracy is the invariance of the equation ( |112| ) under Lorentz 
transformations, under all discrete symmetries of the Dirac equation and, for g = const, 
even scale invariance and translational invariance in momentum space. If, for example, 
0\ = 2 = 75, then, due to the Ward identity 

P^lip, q) = lsG-\ qi ) + G-\q 2 ) l5 , (115) 

= 75G- 1 (g 1 )+G- 1 (g 2 ) 7 5 (116) 



is the divergence of the axial current. Hence, if the equation ( |112| ) is satisfied for p^T 5 ^, it 
has to have the solution (|115|) . 



In order to show that (|112|) has a decreasing solution at \q 2 \ — > oo for any p value, let us 



notice that in Euclidean space this equation has the structure of the Schrodinger equation 

(-d 2 + U)ty=sty (117) 
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in four dimensions at e = 0, with a potential depending on q, spin variables and the external 
vector p. For such a potential the total four- dimensional angular momentum is not conserved, 
only its projection \i onto p v . An equation of this type has always non-singular solutions 
with incoming waves of given /x. E.g., considering an incoming wave with fi = for the 
pseudoscalar 0, we will have a solution 

0o = 75 Co + decreasing scattered waves 

at q 2 — > oo. Or, taking an incoming wave of the form 75^7^ = 75P, we will have 



01 — 7sP + decreasing scattered waves. 

Suppose we found the solution (p 1 . Then, due to the fact that for q — > 00 the solution ( J115|) 
behaves as 

7 5 G- 1 (g 1 ) + G'- 1 (g 2 )75 - Z" 1 + 2 " ^ 



we will find that 



= Z- 1 0i-7 5 G'- 1 (g 1 )-G'- 1 (g 2 )7 5 



2Z- 1 75 z/ 3 



'118 



;ii9) 



is decreasing with | q 2 \ — > 00 at any p. 

By stating that the equation for bound states has a solution at any value of p we do not 
mean that there are no bound states; it means only that the mass of the bound state has to 
be calculated independently. The most natural way to do this is to calculate the self-energy 
of the state 

E(p)= ^xy^v< >wv (120) 



and to solve the equation 



£(p) = 



(121) 



Alternatively, one can calculate the forward Compton scattering of the bound state on a 
fermion, 



(122) 






and after that integrate over the distribution of fermions in the vacuum. But this way we 
will never get a massless Goldstone state. The reason for this is the almost obvious fact that 
( p. 2 Of ) and (|122j) contradict the condition of axial current conservation. 

Let us consider the condition for current conservation in detail. If we introduce f;, as a 
set of diagrams 



f I = 7„75 + 



(123) 
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with a massive fermionic Green's function, it will not satisfy the Ward identity. However, 
the Ward identity will be satisfied by the sum of T 5 and of the Goldstone contribution: 




+ 



12 qi 




(124) 



namely \p = q 1 - q 2 ], 



p M rJ = Pfl fl - ifg = 7B G-\ qi ) + G- 1 (g 2 ) 75 



(125) 



Here g is the Yukawa coupling of a Goldstone boson to a fermion [and / is related to g via 
the fermion loop diagram] 



if Pi* = 7^75 




9- 



(126) 



Knowing P/jT 5 „, we can define the Yukawa coupling g by ( |125| ), ( |126| ). Let us apply the 
operator — d 2 + U to (|125|) . We find that gf satisfies the equation (\112\) since PpT 5 ^ and the 
right-hand-side of ( |125| ) satisfy the same equation. But the existence of ( |124[ ) implies that 
the mass of the Goldstone boson has to be zero. 

In order to clarify the situation with Compton scattering, let us consider the Ward 
identity for the amplitude 



rt(? 2 ,?i,&) 




(127) 



where the dashed line with the momentum k' corresponds to the axial current, the wavy 
line with momentum k corresponds to the Goldstone state, and the other two lines to the 
fermions. The Ward identity for this amplitude is 

k' Si.? 1 k' 



k'^(q 2 , qi ,k) 



75 



g(k,qi) 




+ 




(128) 



g(q2,k) 
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To fulfil ( |128j ), T^(g 2 , qi, k) has to be equal 




+ 




(129) 



A(q 2 ,qi,k',k) 
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here rj^ is defined by ($TZ§. Using (JTgg) and JI2BD we obtain 



i/A = 75^(A;, gi) + #(<?2, fc)7s • 
At large qf, q\ values we will have 

.275 v 3 . Au 3 
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A 



(130) 



(131) 



/ q 2 ' fW 

This means that similarly to the case of asymptotically non-free theories the Goldstone- 
fermion scattering amplitude does not depend on the momentum of the Goldstone boson; it 
decreases only as a function of fermion virtuality. Under these circumstances it is obvious 
that even in an asymptotically free theory the Goldstone boson has kind of a point-like 
structure. 

Amplitudes for the interactions of many Goldstones with fermions can be found in an 
analogous way and have the same property [the property of being independent of the Gold- 
stone momenta]. The self-energy of the Goldstone state is now different from (|120|) . It 
contains two terms 



S(p) 



+ 




(132) 



and equals p 2 at small p 2 . 

As we already have said, in this approach the Ward identity becomes the definition of 
the Yukawa coupling g (the wave function of the Goldstone boson) through p^T 5 ^ which 
has a clear diagrammatic meaning. The equation contains the amplitude / of the transition 
between the Goldstone and the axial current. The expression (|126|) for fp^ is highly symbolic, 
because it contains overlapping divergences. In order to write a sensible expression we can 
use the same procedure as in section ^ where we calculated the polarization operator of the 
photon. Applying the Ward identity, we can write instead of (|126|) 




Differentiating (|133|) over p we will get for p = 



fS^v = 7^75 



(133) 




7^75 



(134) 



[This diagrammatic expression contains exact fermion propagators and bare vertices.] If 
we want to get rid of the 75-s we have to commute 75 with the Green's functions and the 
interaction vertices along one of the fermionic lines. Doing this, we obtain 



{75,G- 1 } 75 




lv + 7m 




{75, G- 1 } {75, G- 1 } 



7m 




(135) 
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Due to the conservation of the vector current, the first two terms in ( |135| ) are zero at p = 0. 
The first one is just the photon polarization operator at p = 0, the second one is the amplitude 
for the decay of a zero momentum scalar into two zero-momentum photons. 

The last term looks like the zero-momentum pseudoscalar-photon scattering amplitude 
which also has to be zero. This, however, is not true, because it does not contain all the 
necessary diagrams. Observing that this term does not contain overlapping divergences we 
can write 

{75,67-!} {75, G- 1 } {75, G- 1 } {75, G" 1 } 

{75, G- 1 } {75, G- 1 } 

In the zeroth order in -, 

f^lJj^^^G-^Gi^G-^GA^A.iq)). (137) 



7 Flavour singlet and flavour non-singlet Goldstone 
states. The U(l) problem 

Up to now we discussed the Goldstone states in a relatively abstract way without fixing the 
concrete asymptotically free theory. In real QCD we have quarks with different flavours and 
there is a difference between flavour singlet and flavour non-singlet states. In order to clarify 
the picture it will be useful to describe the Goldstone state in a different way. 

The previous discussion shows that the Goldstone states in asymptotically free and non- 
free theories are rather similar. Therefore it is natural to try to introduce the Goldstone 
boson in the usual way as a point-like state and to see how this state will interact. In the 
usual discussion of a Goldstone particle we suppose that there is a point-like pseudoscalar 
interaction between this particle and a fermion with the pseudovector coupling 

S =piA- ( 138 ) 

<■ C « Jo 



This interaction induces the radiative correction to the propagator D(p 2 ) of the pseudoscalar 

D(p 2 ) = 'vvvw^ + 'v<^^>\^ + ... , (139) 

where the wavy line describes the bare pseudoscalar Green's function Dq. If the fermion is 
massless and the axial current is conserved, this pseudoscalar will not interact; its self-energy 




s b ) = 7f?W75< PPuT u = PiJhtTp, (140) 
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is equal to zero. If, due to symmetry breaking, the fermion becomes massive, it starts to 
interact and acquires a self-energy different from zero. By using the diagrammatic definition 



of It [see (p3Dl, we will find 



2 f 2 



11 

/o 2 



(141) 



where / is the same amplitude for the Goldstone-current transition as what was discussed 
in the previous section. Hence, 



D(p 2 



Do 1 f$ - P 2 f 



In the limit /o — > we will have 



D{p 2 



f 2 
J o 



P 2 f 2 



(142) 



(143) 



The pseudovector interaction between Goldstone and fermion is jP^r^ with a pseudovector 
coupling j defined by the fermion mass. The limiting procedure /q — > can be understood 
if we accept that the interaction responsible for symmetry breaking changes the fermion 
vacuum fluctuations not only at finite momenta but also near the ultraviolet cutoff. This 
change in the fermion vacuum fluctuations is responsible for the leakage of the axial current 
in [from] the region of finite momenta; it can produce the driving term for the Goldstone 
state, recovering the current conservation. 

In general, the pseudovector coupling has a disadvantage compared to the pseudoscalar 
coupling which we have discussed before: it looks unrenormalizable. But in the case of 
flavour non-singlet states it can always be replaced by the pseudoscalar coupling with the 
help of the trivial relation [p = qi — q?\ 



Pl5 = (<?2 - ™(<?2))75 + 75(?i - + ^(92)75 + l&n(qi] 



(144) 



which leads to the Ward identity ( |125f) for pseudoscalar coupling. If we include the emission 
and the absorption of Goldstone bosons inside the diagram, then in the process of this 
replacement a point-like amplitude appears which corresponds to the quark interaction with 
many Goldstone bosons. Nevertheless I believe it is possible to prove that, due to the 
decrease of these amplitude as functions of quark virtuality, the theory is renormalizable. 

Due to the anomaly, the case of a flavour singlet current is very different even on the level 
of the Goldstone Green's function. In this case the corresponding polarization operator E 
will contain not only the quark loop which we have discussed but also a gluonic contribution 



S(p 2 





(145) 



The triangle diagram f^ v present in ( 145|) was calculated many years ago by Adler, Bell and 
Jackiw [3] 



a 



7T 



(146) 
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With this expression for f^ V) S(p 2 ) still has the form ( p. 41 ) but it will be quadratically 
divergent and 



a 



/«+("] A 



2 

UV 



oo 



(147) 



where Auv is the ultraviolet cutoff. This means that Goldstone particles exist in the anoma- 
lous case but they are decoupled from all physical states. At the same time the Ward identity 
still makes sense because the product gf does not depend on /. Nevertheless, the concrete 
form of the Ward identity will change. The reason for this is again the Adler-Bell-Jackiw 
anomaly [3]. For the triangle diagram 



per 




(148) 



the replacement of the pseudovector coupling by the pseudoscalar coupling gives an incorrect 
result: instead of the correct expression 



pa 




a 



71 



(149) 



which was obtained in [3], we get just the first term. We can try to write the Ward identity 
using ( |149| ), but this turns out not to be necessary. The reason is that in this approach the 
axial current between gluonic states, defined symbolically by the relation 



■PP _ "TP i 

per pa 




(150) 



(where the term is defined diagrammatically and g pa = ^jY v pa ) is just the transverse part 

V'L- (151) 



of 



r p = y 

per per 



I'pl'l' pi/ 



2 Per' 



The axial current between quark states can be written as 



~p 



r^ = r (g 2 ,gi)+i% 2 , gi ). 



(152) 



Here 



+ 



(153) 
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is the same set of diagrams as in the non-singlet case and is the " axial current of gluons' : 



9 



rp 

6 

^ — ^ — 



+ 



ip 

6 
1 1 \ 
/ 1 \ 



+ . . . 



(154) 
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In the same way the Goldstone-quark interaction can also be divided into two parts. In 
the first part we can replace the pseudovector coupling by the pseudoscalar coupling. The 
second part is the longitudinal part of T^. Consequently, 




9 ' 



(155) 



and the Ward identity can be written in the form 



p,T -ifg = l5 G-\q 1 ) + G- 1 {q 2 ) lb . 



(156) 



This expression enables us to answer the question, what happens with particles like rj' which 
would be Goldstone states if we did not take into account that they can decay into two 
gluons. Let us consider the contribution of the massive pseudoscalar flavour singlet particle 

r]' to the Ward identity ( |156| ). The divergence p^JP has a pole corresponding to r]': 




ip 2 f, 



V 2 2 



fly 



(157) 



The Yukawa coupling of the Goldstone state has also a pole: 

{75, G- 1 } 

-'7.7= >' /!/ ' {':,(■■ 



/r — P A 



158) 



sy 



The right-hand side of ( |156| ) has, however, no poles. This condition can be satisfied if 

l?U = {^75, G-\q)} <^ 9rf . (159) 
The same Ward identity ( |156| ) [substituting ( |157| ), ( |158| ) and using ( |159| )1 gives at p 2 = 

frfQif = {ils,G-\q)} , (160) 
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from which equation it follows that \j? for 7/ is equalQ 



(161) 



This means that rf acquired a mass due to the transition into a Goldstone boson which itself 
is decoupled. 



It is interesting to notice that at relatively small fi 2 when the comparison between Ward 
identities for different values of p 2 (p 2 = fi 2 , p 2 = 0) makes sense, ( |161| ) gives us the same 
result as the equation (|120|) (without the additional point-like term which is present in the 
7r-meson case ( |132| )). Indeed, in ( |120| ) we can write 






0. 



The sum of the first two terms is equal to p 2 (as in ( |132j )) and therefore 



p =/i 




(162) 



p=0 



in agreement with ( |161| )ffi We see that the subtraction term in the 7r-meson self-energy, 
reflecting the quasi-local structure of the pion, disappears in the case of rj' '. In this sense 77' 
is a normal bound state without a point-like structure. 

The approach we presented here for the resolution of the [/(l)-problem is technically 
very close to the approach developed by Veneziano [4] but the underlying physics differs 
essentially. In Veneziano's approach big long-range fluctuations are responsible for the 7/ 
mass. In our approach 77' is a normal qq bound state with no local structure which would be 
responsible for its small mass if it were a Goldstone state. This local structure is destroyed 
by the decay on hard gluons. 



8 QCD with massive quarks 

We have discussed in detail an asymptotically free theory with massless fermions. We came 
to the conclusion that in order to obtain the correct spectrum of Goldstone particles, axial 
current conservation has to be imposed on the theory. QCD, however, contains massive 
quarks and the same spectrum of massive quasi-Goldstone particles (the pseudo-scalar octet) 
as the theory with massless quarks. The problem is, how to impose the condition of axial 
current conservation on a theory which obviously does not conserve the axial current. 

Strictly speaking, I don't know how to do this. For our real world, however, there is a 
natural possibility to solve the problem. 

3 [The overall sign in the following equation differs from that of the original manuscript.] 

4 [Notice however that the expressions (161) and (162) have different signs; see the previous footnote.] 
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In the real world QCD is part of the standard model describing strong, electromagnetic 
and weak interactions. In the standard model all particles are supposed to be intrinsically 
massless and their masses appear as the result of symmetry breaking due to some kind of 
Higgs mechanism with or without elementary Higgs particles. The possibility of such a 
mechanism is guaranteed by the conservation of the left-handed S"?7(2)-current j®. For the 
matrix elements of this current between any two quarks with momenta q 2 , qi, we have 
the Ward identity [p = q\ — q 2 ] 



1 + 75 T a 



G-\qx) 



(163) 



In the case of massive fermions p^T^ contains the contribution of three Goldstone bosons 
responsible for the masses of W ± and Z°: 



p^ a =p^ a -fg\ [9 a = j9\- (164) 

For large g 2 -values G~ l contains a massive term Z~ 1 mQ. Hence, at large q 2 we have 

(f9 a )o = ^{Ta^moZ- 1 } + ^moZ- 1 ]. (165) 

At small q 2 of the order of the QCD scale, fg is defined by the total quark mass m(q) which 
for the light quark is much larger than m . Using the relations ( |163j ) - ( |165| ) we can calculate 
the masses of the pseudoscalar octet in the same way as we did for r]'. 

Suppose there is a bound state of light quarks which is a pseudoscalar particle (the 7r- 
meson) with a finite mass /i. The pole corresponding to this particle will contribute to both 
terms in (164). With this pole contribution, p^T^ and fg can be written in the form [with 
np standing for "non-pion" or "non-pole"] 




UP 2 - 2 € + (P^aU, 

/r — p z 



(166) 



fg a 



(fg)o 




+ (fg a ) np = (fg a W 



>9l , 



^ — p* 



A + Ug a \ 



np ■ 



(167) 



As in the case of r]', the condition for the pole cancellation gives us the value of /i : 



2 x 

/i • o ab 



T (fg a )o< 

Jn 



9v 



(168) 
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At p 2 = we have 

Ugl + (fg a ) np = - A {r al5 ,m(q)Z- 1 } + \ [r^mZ- 1 ] . (169) 

[Here m is the light-quark mass matrix, m = m s + \[m u — m^r^, with m s = \(m u + ma)-] 
The Yukawa coupling g n of the bound state has to decrease at large q 2 . Therefore we have 
to identify (fg) np with (fg) . Then 

Ug% = ^j 5 m s (q)Z-\ (170) 

where m s (q) is the isotopically invariant part of the quark mass term behaving at q —>■ oo 
as v 3 jq 2 . Finally, (fg)o oc mo is defined by weak interactions with Goldstone states, and 
f^g^ oc m s is determined by pion exchange. The conclusion of these considerations is that 
in the case of massive quarks the conservation of the left-handed SU(2) current can play 
the same role for the calculation of the coupling and the masses of the flavour non-singlet 
pseudoscalar particles as in the massless case the conservation of the axial current. The 
result is that the masses become different from zero while the Yukawa coupling remains the 
same at least for small mo values. 

The expression ( |168| ) for fi 2 can be also obtained without discussing the current conser- 
vation, just by calculating as it is defined by (|132|) for the quark mass m s + m , keeping 
only the term linear in m . It is in agreement with the idea that m affects the mass of the 
Goldstone state but not its wave function, which always has been the common belief. 

It is not trivial that, if we calculate /i 2 by using ( |168| ) or ( |132| ), we get a logarithmic 
divergence. For the 7r-meson mass the main divergent part equals 

m '^r ; LT-f mMAq) (m) 

(here v Q is of the order of Aqcd)- 

In the standard model the behaviour of mo and v 3 can be calculated in a fantastically 
wide region of q 2 : from 1 GeV up to the scale where one of the couplings of the standard 
model (the ?7(l)-coupling gi, the Yukawa coupling h or the coupling of the self-interaction 
of the Higgs particles, A) becomes of the order of unity. If the [/(l)-coupling gi is the first 
to become large, which seems natural, this scale is Auv — 10 38 GeV. 

At q 2 > A^Y the behaviour of mo and v 3 is unknown. Because of this, m n is not calculable 
in principle and it has to be considered as an arbitrary parameter. If we, however, assume 
that at q 2 larger than Auv, m o(l) and v 3 (q) vanish, we will be able to calculate m^ and, 
surprisingly, this calculation gives a reasonable value for m n with Auv — 10 38 GeV [2]. The 
expression (|171|) for m 2 is also in agreement with the naive expression 

9 2m n , - . 
ml = 

J 7T 

The important difference is that now (^^/) is determined not only by strong but also by 
weak interactions: it goes to infinity if the weak interaction is removed. 
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9 The pion contribution to the equation for light-quark 
Green's functions 



From the previous discussion it is obvious that the small-mass pion contribution has to be 
included in the equation for the light-quark Green's function. Fortunately this is very easy to 
do. Having in mind that now the diagrams contain not only the gluon contribution but also 
the emission and the absorption of pions we will find that as before, the main contribution 
to d 2 G~ 1 comes from the simplest diagram r J~ yv/v \ with the coupling {275, G -1 } at zero 
momentum instead of the gluon coupling <9 M G _1 . It leads to the following equation for the 
Green's function: 

&G~\q) = gd,G-\q)G(q)d,G-\q) - { ll ^ G ~ l }G(q){l l5l G- 1 } T7 ^ (172) 



The equation for bound states 
diagrams 



has also to be changed. The correction comes from the 



7T 






Instead of (1271) we will have 



d 2 r(p,q) 



+ 



g{q){A v {q 2 )d v <t> a (p, q) + d„<p a (p, q)A v { qi ) - A,(g 2 )0>, q)A l/ {q 1 )} 
~~7~2~~c2 h^G-^jGiq^r^Giq^i^G-^)} - 



XG(q 1 ){i lB ,G- 1 (q 1 )}^- - { tl5 , G ^(q 2 )}G(q 2 )^\ 



(173) 



Here A is the emission amplitude of the zero momentum pion in the transition of the bound 
state to the gg-pair. This amplitude has to be defined by the axial current conservation. 
There is another important quantity in this equation, namely: /„-. In section |6| we have 
written an explicit expression for f 2 including only the gauge field contribution and ignoring 
the pion contribution. Now we include the pion contribution in the equation for the Green's 
function and, to be self-consistent, we have to do the same for f 2 . I was not able to carry 
out this in any order in the Yukawa coupling, but in the first order in g 2 the equation ( |136| ) 
is correct if one adds the diagrams 



d^G- 1 



d^G- 1 



{75, G- 1 }<^W^0> {75, G- 1 } + {75, G- x }< < > {75, G- 1 }. 




(174) 



As we see, the gluonic correction of the order of ^ and the pionic correction of the order of g 2 
have the same diagrammatic structure. In order to estimate the value of these corrections, let 
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us take just the contributions of zero momentum gluons and pions to them. It can be shown 
that the contribution of zero momentum gluons cancels in the last two diagrams of ( |136| ). 
Zero momentum pion contribution comes only from the first diagram of (|174|) . Transferring 
the differentiation from the fermionic line to the pionic line in this diagram we obtain the 
contribution of zero momentum pions in the form 

{75, G~ 1 }<^7^J^> { 75 , G- 1 }^ . (175) 

The expression for which includes the zero-momentum pion contribution is 

Sf, = j-^ ri T<:({n,,G- 1 }G{i 7 ,,G- 1 }GA ll A^ 

+ shlWfi^^' ^- (176) 

It gives us an understanding of the scale of possible pion contributions. It is interesting to 
note, that (|176|) is not an expression in terms of the Green's functions but an equation for 

In the next paper we will show that the pion contribution changes essentially the structure 
of the equation for the Green's function. The new equation has a solution corresponding to 
the confined quark. At the same time the symmetry breaking solution will not necessarily 
survive (at least if g(0) is large). 
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